ICRR-Report-566 



WIMP dark matter in gauge-mediated SUSY breaking models 

and its phenomenology 



Junji Hisano 0- ' 6 , Kazunori Nakayama", Shohei Sugiyama a , Tomohiro Takesako", 

and Masato Yamanaka a 

a Institute for Cosmic Ray Research, University of Tokyo, Kashiwa 277-8582, Japan 
b Institute for the Physics and Mathematics of the Universe, University of Tokyo, 

Kashiwa 277-8568, Japan 

Abstract 

We propose an extended version of the gauge-mediated SUSY breaking models 
where extra SU(2)l doublets and singlet field are introduced. These fields are 
assumed to be parity-odd under an additional matter parity. In this model, the 
lightest parity-odd particle among them would be dark matter in the Universe. In 
this paper, we discuss direct detection of the dark matter and the collider signatures 
of the model. 



1 Introduction 



Supersymmetry (SUSY) is one of the best-motivated physics beyond the standard model 
(SM) since it naturally solves the gauge-hierarchy problem. SUSY must be spontaneously 
broken, however, in order for SUSY particles to obtain sizable masses. Among proposed 
SUSY breaking models, the gauge-mediated SUSY breaking (GMSB) models pQ are in- 
teresting, since the SUSY flavor problem does not arise in GMSB models [2]. 

Although low-energy phenomenology of GMSB models seems to be quite successful, it 
has non-trivial aspects from a cosmological point of view. In GMSB models the lightest 
SUSY particle (LSP) is the gravitino, superpartner of the graviton. The relic abundance of 
the gravitino in the Universe depends on the reheating temperature after inflation, and a 
stringent upper bound on the reheating temperature is obtained so that the gravitino does 
not exceed the present dark matter (DM) abundance [3]. In particular, this implies that 
GMSB models are not compatible with the thermal leptogenesis scenario [I] for most 
range of the gravitino mass. An exception is the low-energy GMSB models where the 
gravitino is lighter than about 10 eV and no upper bound on the reheating temperature is 
imposed [5J E]. In this case, no candidate for DM exists in the minimal setup of GMSB0 

One may consider that the QCD axion [7] can play a role of DM. However, cosmology 
of SUSY axion models is quite non-trivial taking into account the existence of axino, 
fermionic superpartner of the axion, and saxion, scalar partner of the axion. The axino 
is produced thermally in the early Universe with a significant amount, and hence the 
reheating temperature is more severely restricted [HI El HO]- The saxion coherent oscillation 
and its subsequent decay also gives catastrophic cosmological effects unless the reheating 
temperature is sufficiently low [HI [12]. Thus to make the axion the dominant component 
of DM in GMSB model requires careful considerations^] 

In this paper we extend the GMSB models to include a candidate for WIMP DM. 
The minimum extension might be to add chiral supermultiplets with fundamental repre- 
sentation of SU(2)l (H', H'). Those fields are parity-odd under an additional Z 2 -parity 
assigned. If all other minimal SUSY standard model (MSSM) fields are even under the 
.^-parity, H' and H' can be stable and a candidate for WIMP DM. Since they are weakly- 
interacting, their relic abundance falls into a correct range in the thermal history of the 
Universe. However, this kind of extension is already excluded since it predicts too large 
direct detection rates in the DM search experiments through the coherent vector coupling 
to nucleons by Z-boson exchange. In order to avoid the direct detection bounds, we fur- 
ther introduce a singlet chiral multiplet S' which is also parity-odd. Then WIMP DM 
becomes a mixture of S', H' and H', and the lightest parity-odd particle is real scalar bo- 
son or Majorana fermion so that the vector coupling of DM-nucleon is forbidden. Sizable 

1 In the low-energy GMSB scenario, a model in which a baryonic bound state of strongly interacting 
messenger particles becomes cold DM is proposed [13l[T4]. In this case DM has a mass of 0(100) TeV, 
and cannot be detected by accelerator searches and direct detection experiments. 

2 The strong CP problem may be solved in the low-energy GMSB models through the Nelson-Barr 
mechanism [TSl US] . When the SUSY breaking sector is decoupled with the spontaneous CP-violating 
sector, the radiative correction to the QCD-thcta term is suppressed due to the non-renormalization 
theorem. Thus the axion is not necessarily needed in this case. 
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Table 1: Particle contents of the model. 



SU(3) C SU(2) L U(1) Y Z 2 



H 
H 
H' 
H' 
K 

K 



s' 



1 2 -1/2 even 

1 2 +1/2 even 

1 2 -1/2 odd 

1 2 +1/2 odd 

3 1 +1/3 odd 

3 1 -1/3 odd 

1 1 odd 



interactions still exist through the Higgs exchange process, which is within the reach of 
on-going/future direct detection experiments. 

This paper is organized as follows. In Sec. [2] we define our model and study the 
properties of the DM particle. In Sec. [3] the DM direct detection rate is evaluated and it 
is shown that it is within the reach of current/future direct detection experiments. We 
discuss the LHC signature of this setup in Sec. HI Finally in Sec. |5]we give conclusions. 



In this section we define our model and discuss the properties of the DM particles, such 
as their mass, spin and interactions, in the model. We show parity-odd chiral multiplets 
newly introduced and MSSM Higgs doublets in Table [TJ The lightest particle among 
mixtures of parity-odd fields is stable and can be a WIMP DM candidate. In order to 
maintain the unification of the gauge couplings, we also introduce SU(3)c triplets (H' c , 
H' c ) which compose 5 and 5 of 377(5) with H' and H'. The most general renormalizable 
superpotential is 



SU(2)l products are defined as HH = H°H° — H~H + . Five real parameters and one 
CP violating phase, 9 = arg(/i// M s A*A£) , are introduced. For simplicity, we take 9 = 0, 
so that all parameters including /i' are real and positveJl 

3 The electric dipole moments (EDMs) are induced by electroweak two-loop diagrams in our model. In 
Ref. [17] the EDMs are discussed in a model similar to ours, where SU(2)l doublet and singlet fermions 
are introduced in the standard model. It is found from their result that when the CP violating phase and 
the couplings are 0(1), the electron EDM induced by parity-odd fermions reaches current experimental 
bound for fj,' ~ M s IS 1 TeV. When M s < //, the constraints are more loosened. Thus, the constraints 
are not severe at present, though future EDM searches might give severer constraints on the model or 
find the signature. 
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The n parameter in Eq. (P) is of the electroweak scale from the naturalness argument. 
In addition, the mass parameters fjf and M s are also expected to be of the electroweak 
scale so that the DM relic abundance explains the WMAP measurement. Those mass 
parameters should have a common origin. In the low-energy GMSB models, an extra 
dynamical sector is introduced to generate a constant term in the superpotential for the 
cosmological constant to vanish. This sector can also produce the dimensional couplings 
in Eq. (JTJ with the magnitude of 0(100) GeV [Tg]0 

Soft SUSY-breaking terms for the parity-odd fields are generically given as 

Vsoft = m 2 H ,\H'\ 2 + m\,\H'\ 2 + m 2 s ,\S'\ 2 

+ (B'^H'H' + ^B' S S' 2 + A^H'HS' + A^HH'S' + h.c). (2) 



We study particle spectrum based on the minimal gauge-mediation model (MGM) 
throughout this work, in which the messenger sector is composed of vector-like 5+5 
representations of 577(5). In the MGM, A- and 5-terms vanish at the messenger scale 
M mess . In this model, the scalar squared mass of the singlet is also zero at M mess . The 
only non-vanishing terms at M mess in Eq. (jSJ) are 



m H' — m H' — 2iV mess 



1 ray\ 2 3 /«2\ 2 ' 
4V4vJ 4 UtJ 




where N meas represents the number of S77(5) representations introduced as messenger 
fields, A = (F)/M mess , and (F) is the vacuum expectation value of the F-term which 
couples to messenger fields. The loop function f(x) is given in Ref. [2]. 

Now we discuss the properties of DM particles. The mass matrix for parity-odd 
fermions in this model is 

-£=^ T M F i> + h.c, 7/>=|/7' u |, (4) 



-Ms —Xiv — \2v\ 
Mp = | -Ait; -fi' , (5) 
-\ 2 v -/I 1 / 

where v = (H°), v = (H°), and v 2 + v 2 = 2m 2 z /(g 2 + g' 2 ). Similarly, the squared mass 
matrices for parity-odd bosonic states are 

-C= \vlM { : )2 VR + ^JM«Vi, < PR =\H$\,<p I =\H?\, (6) 







fs'A 


K 


1,^=1 











4 Fortunately, the messengers do not generate the S-terms in this mechanism at the leading order 
as mentioned in text, and we do not need to care about the so-called /i/-B/i problem. When the Higgs 
multiplets have (direct or indirect) couplings to SUSY-breaking field to generate the /i term, the term 
would be too large. 
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Figure 1: Mass, spin and type of the DM particle on n'-M$ plane. Dotted blue lines are 
contours of DM mass. DM is boson below the solid red line. It is fermion above the line. 
DM is doublet-like on the left-hand side of the green broken line, and it is singlet-like 
on the right-hand side of the line. The input parameters are given in text. Gray bands 
are regions where the relic DM abundance is consistent with the WMAP result within 2a 
level. 



tX\v 2 + X\v 2 + Af| + m 2 s , ± B' s X 2 fJ.'v + XiMgv ± (X 1 fiv + A*v) X^'v + A 2 A/ S t; ± (X 2 fiv + A$v)\ 
M^ )2 = lx 2 n'v + X 1 M s v±{Xinv + A^v) /t' 2 + X 2 v 2 + m\, + A X 1 X 2 vv ± 

\Xin'v + X 2 M s v± {X2HV + A^v) XiX 2 vv ± B'^ fi' 2 + X\v 2 + m\, + A / 

(7) 

A's are contributions of D-term potential; 

A = ^m|cos2/3, A = -A, (8) 

where tan/3 = v/v. The fields with subscripts R and / are CP-even and odd states, 
respectively, and canonically normalized as = -^(0_r + i(j>i)- 

We show mass, spin and type of the lightest parity-odd particle in Fig. [TJ In the 
numerical calculation, we assume MGM with 5 + 5 messengers, and we use the result of 
[20] . Here we took Ai = A2 = 0.3, tan/3 = 42, \i = 660 GeV, N mess = 1, the gluino mass 
1 TeV, and A/M mess = 0.5. fi and tan/3 are fixed from electroweak symmetry breaking 
conditions and a condition that = at the messenger scale. We set these value as 
reference point throughout this work. 

If jjl 3> M$, DM is singlet-like CP-even boson. The singlet bosons get no SUSY 
breaking mass terms in the GMSB model at the leading order. The bosonic and fermionic 
states are degenerate in masses. In our set up, the F-components of H and H generate 
mass splitting between CP-even and odd states. Thus, one of the bosonic states tends to 
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be lighter than fermionic one. (If we take y! negative and |//| 3> Ms, DM is singlet-like 
CP-odd boson. In this case, mass splitting source of CP-odd boson is larger than that 
of CP-even boson as expected from Eq. ([7]).) If Ms is too small, the lightest bosonic 
state becomes tachyonic and the ^-parity is broken spontaneously. If // ~ Ms, DM is 
a fermion which is mixture of singlet and doublets. In the region of /i' <C Ms, DM is 
doublet-like fermion. The doublet-like bosonic states are heavier than fermionic one due 
to the GMSB effect. 

Now we see that this model predicts a correct relic DM abundance measured by 
WMAP [2TJ with appropriate parameters. Since our model includes particles with mass 
close to the DM mass, we should take into account the coannihilation effect for calculating 
the relic DM abundance. 

First, consider the case where DM is singlet-like CP-even boson (S' R ). In this case 
we consider three coannihilating states {S' R , S'j, S'). The main annihilation processes in a 
non-relativistic limit are S' R S' R — > h°h° and S'jS'j — > h°h°, where h° stands for the SM-like 
Higgs boson, as far as they are kinematically allowed. Here we can neglect other processes 
like S'S' h°h°, S' R S' ->• h°H°, and S' R S'j h°A° where A is CP-odd Higgs boson, 
and so on, because H° and A are assumed to be heavy. The reason for omitting the S' 
annihilation is as follows. Since S' is a Majorana fermion, in the non-relativistic limit, 
a pair of them is in CP-odd state. Then the opening channel from the viewpoint of CP 
conservation is the pair annihilation of S' into a pair of h° in p-wave state. This channel, 
however, is forbidden by the angular momentum conservation. 

The DM effective cross section [22] of these processes is given by 



where 



= ( TT^) <™>s^->ft°w> + (j^ ) M^-ww, (9) 



1 1 fm h0 y 2 (M 5 Ax + ^) 2 Gu'Ax + z^ 



„£>„^ri"n" r- a - 2 \ I /I I -2 i -2 2 -2 -2 2 

bArcrrig, y \ m s' / \ 171 h' + m s' ~ m h° m H i-\-m A s , — m^ 

(10)' 

(M 5 Ax-^) 2 (//% - //A 2 ) 2 




647r m 2 ^ \] \ m s'j J V 1 + - m lo m 2 H ^ + m |j - m 2 h o 



(11) 

_ 3 _ _ 

w = ( — ^- ) exp -x— x - ^ , (12) 

\ m S' R J V / 

x = 7 -^^. (13) 

freeze out 

Here, ms' R , m s'j, r nH[, m H' 2 are the mass eigenvalues and Tf reez e out is the freeze-out tem- 
perature of S' R . The relic DM abundance, in terms of the density parameter Qs' R h 2 , is 



expressed as 



,'2.7 x 1(T 9 GeV- 2 \ 
n s , R h 2 = 0.10 . . 14 
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The gray region in Fig. 1 (bottom one) shows the parameter space in which the DM 
relic abundance is consistent with WMAP observations within 2a level [21] ■ Thus we 
can see that the relic abundance falls into the correct range measured by WMAP for 
certain choice of parameters. We note that the coannihilation effect can change the relic 
abundance only 10% or so in the correct range for DM abundance.® 

In the case of doublet-like fermion DM, we consider four coannihilating states ( 

H'° , H'°, H'~ , H' + ). The cross section and its relic abundance are basically the same 
as those of MSSM Higgsino-like DM whose mass is heavier than VT-boson mass. If the 
coannihilation effect is efficient, the DM effective cross section and its relic abundance are 
given in Ref. [23] ; 

o 4 

(a cS v) = — f - (21 + 3 tan 2 6 W + 11 tan 4 6 W ) , (15) 



5127T/X- 



2 



n A , oh 2 = o.io (j^J , (16) 

where four states are taken to be degenerate in massif] The effective cross section in 
Eq. ( 1T5|) drops final states including ordinary MSSM superpartners, because their total 
contribution is at most 10% in our reference point comparing with the SM final state 
contribution. The correct DM relic abundance consistent with WMAP observations within 
2cr level is obtained in the gray region in Fig. 1 (top one). 

Finally, we discuss about H' c and H' c in Eq. (OQ), which are introduced in order to 
maintain the successful gauge coupling unification. They have color and electric charges, 
and their masses are also expected to be order of electroweak scaled In the model, they 
are stable at the renormalizable level and unstable due to the higher dimensional operators 
suppressed by relevant physics scale M p h ys such as 

tt^(55sm)(55sm). (17) 

^ W phys 

In the GMSB model, triplet fermion is lighter than bosonic one. When DM is doublet- 
like fermion, the lifetime is 0(1) — O(10 3 ) sec for masses 100 — 1000 GeV in the case 
that M phys = M G ut — 10 16 GeV. Here, it is assumed that the triplet fermion can decay 
into the MSSM SUSY particles directly. When DM is singlet-like boson, the decay of 
the triplet fermion is suppressed by the mixing between the singlet and doublet states. 
The feature of these long-lived colored particles are constrained from the viewpoint of 
cosmology. Their relic abundance would be determined by geometrical cross sections of 
order of 10 mb effectively at temperatures below the QCD phase transition, and they an- 
nihilate efficiently before the Big Bang nucleosynthesis (BBN) if their masses are lighter 
than about 10 11 GeV EU. Their relic abundance after this annihilation is estimated as 



5 The it final state via stop exchange is ignored in Ref. [23 , and our case does not include such a 
process. 

6 If [i! c is equal to // at the GUT scale, \J C is about twice larger than [i! at the electroweak scale. 
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KTVlX/TeV) 1 / 2 , where rib is the baryon number densityQ The future collider experi- 
ments may prove existence of the colored particles. We discuss its possibility in Sec. SJ 

3 Direct detection rates 

In this section we evaluate the DM-nucleon scattering cross section in our model. DM 
particles in the model have sizable interactions with nucleons and hence they may be 
detected through on-going or future direct detection experiments. 

First let us consider the case where the DM particle is fermionic (denoted by x). DM 
interacts with nucleons through the Z-boson and Higgs-boson exchange diagrams. The 
former yields a spin-dependent (SD) and the latter yields a spin- independent (SI) effective 
interactions. The effective Lagrangian for these interactions is written at the parton level 
as 

£ e ff = Yl d iXYlhXQl^im + fiXXQQ, (18) 

q=u,d,s q=u,d,s,c,b,t 

where coupling constants are given by 

d « = St ((0f)L - { ° F)li) ' (i9) 

and 

q 9 2m w \ m 2 h0 m 2 H0 J ' 

with 

cos a sin a 

C-huu : 77) CHuu '. 77' \ ) 

sin p sin p 

sin a cos a . . 

Chdd = 77, CHdd = 77, (22) 

cos p cos p 

Chxx = ( A i (°f) ii (0 F ) 12 cos a - X 2 (Of) 11 (0 F ) n sin a) , (23) 
1 

chxx = (Ai(0 F )n(0 F )i 2 sina + A 2 (O f ) u (O f ) 13 cos a) . (24) 

Here F is a 3 x 3 mass diagonalizing matrix, which is obtained from Eq. (JSJ). 

Using these couplings, the SI scattering cross section between DM and nucleus with 
mass rriT is expressed as [261 EZj 

A / \ 2 



gffl = - T K/p + ^n/n) 2 , (25) 



7 It is also argued in Ref. [3S] that even such a small abundance may affect the BBN and a lifetime of 
the colored particles should be shorter than about 200 sec. This would give a lower bound on the triplet 
fcrmion mass, depending on -M p hy S and the main decay mode, though we do not include this constraint 
in this paper. 
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and the SD cross section is given by 



0"SD 



7T Krriy + rriT 



T \ P 

(S p ) + a n (S n )) 



(26) 



Here n p (n n ) is the number of proton (neutron) in the target nucleus. J is the total nuclear 
spin, and (S p / n )) = (A\S p r n \\A) are the expectation values of the spin content of the proton 
and neutron groups within the nucleus A [28J, and 



a N = d i Aq 

q=u,d,s 



N i 



2 Stl Aq N = (N\qj^ 5 q\N). 



(27) 



(28) 



where s M is the nucleon's spin and N = n,p. The DM-nucleon effective coupling is 
constructed from the DM-quark effective coupling as follows [29J, 



N 



E 



f AW 
JqJT a 



m N ' m q 

q=u,d,s * 



+ 4/t g Yl It 



m„ 

q=c,b,t 



f: 



...=i- e /r 

?(p) 



(29) 
(30) 



For the nucleon mass matrix elements, we take f^J = 0.023, f™ = 0.034, /£° = 0.019, 

/W = 0.041 [301 |3T] and = /£ } = 0.025 [32]. 

Next let us consider the case of bosonic DM (%). The effective Lagrangian through 
the Higgs boson exchange diagram is written as 



q=u,d,s,c,b,t 

This yields the following SI scattering cross section, 



(31) 



0"SI 



m x m T \ , , ^ 



(32) 



where f p (f n ) is given by Eq. ff29|) . with c/j XX (ch- xx ) lY1 fg (Eq- (l20|) ) is replaced by the 
following, 

c hxx = ((0 B YiO|)n sin a - (0 s Y 2 Ob)h cos a) , (33) 



ch 



XX 



2V2 
1 

2V2 



(-{0 B Y 1 T B ) ll cosa - (0 B Y 2 Ol) n sma) 











(34) 
(35) 
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Figure 2: Contours of SI (top) and SD (bottom) cross sections with a proton on fjf-Ms 
plane. Broken lines are boundaries between fermionic and bosonic DM. 
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Here Ob is a 3 x 3 mass diagonalizing matrix for the mass matrix M B (= M ( B +) ), which 
is given in Eq. ([7]). 

Fig. [2] shows contours of the SI and SD interactions with a proton on fi'-Ms plane. 
Parameters are set to be the same as those in Fig. [TJ The most stringent bound on 
the SI cross section comes from CDMS-II results [33], [33]. The bound reads asi/m x < 

3 x 10" 46 cm 2 /GeV for m x > 100 GeV. The observed two DM-like events at CDMS- 
II [3l] can be explained for appropriate parameters, if we take the two events seriously 
and assume they are caused by DM-nucleon scatterings. The predicted SI cross section 
is within the reach of future or on-going direct detection experiments. 

Among the direct detection experiments, the best bound on the SD cross section 
comes from the XENON experiment |35j (osd/t^x ^ 3 x 10~ 40 cm 2 /GeV for m x > 
100 GeV). In addition, the SD cross section is also constrained by observations of energetic 
neutrinos from the Sun produced by annihilations of DM particles captured by the Sun 
[351 EU ESI W\- Super-Kamiokande @Q], AMANDA [UJ and IceCube with 22 strings 
give the stringent limits [12]. The fermionic DM in the present model mainly annihilates 
into W + W~ , and their subsequent decay modes W — y ev produce high-energy neutrinos 
which may be detectable at IceCube. The current IceCube bound for mass at 250 GeV 
is osd ~ 3 x 10~ 40 cm 2 , assuming the W + W~ final state. 

The bound on the SD cross section is still far from the prediction. However, the 
experiments sensitive to it would be important to determine spin of the DM particle after 
the DM is discovered. Furthermore, we might reconstruct the model by using results of 
the direct DM searches sensitive to SD and SI cross sections under assumption of the 
thermal relic scenario, when the DM is fermionic. This is because the model has only 
four input parameters, //, M s and Ai/ 2 . 

4 Collider signatures 

One of the handles to confirm this model is discovery for the signatures of the extra 
particles at the LHC. In the confirmation, the key issue is the selection of missing energy 
events of this model under the backgrounds of ordinary GMSB models. In this section, we 
estimate the discovery reach of extra parity-odd particles, and discuss the feasibilities for 
the selection assuming that the GMSB model has been already experimentally established 
as the SUSY breaking model. 

First, we discuss the case that the lightest parity-odd particle is doublet-like neutral 
fermion, H'®. The mass difference between H'® and heavier doublet-like neutral (charged) 
fermion, H B (H'^), is less than 40 GeV when \J > 94 GeV within the parameters of 
Fig. [TJ Accordingly, the two-body decays of H'^ and i7 /=t are inaccessible. 

Clean multilepton events, pp — > H'^H !± — > + I'vtiH'®, offer a promising way for 
the selection. In the low-energy GMSB models the SUSY events accompany energetic 
tau leptons or photons since the next-lightest SUSY particle is typically stau or Bino-like 
neutralino. Thus, an observation of the missing energy events accompanying trileptons 
clearly point to H^H'^ production as its source and would be the evidence of underlying 
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Figure 3: Cross sections for pp — > H'^H^ and pp — > H'^H^ . Here H'° and H r± stand for 
neutral CP-even (CP-odd) and charged extra boson, which are components of common extra 
doublet. H'jj and H'^ stand for superpartners of neutral and charged extra doublet boson. We 
took only direct pair production as their production processes, and included all combinations of 
final state pair in each result. 
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Figure 4: Production cross sections for extra triplet boson H' c and its superpartner H' c as a 
function of their mass. We took only direct pair production as their production processes. 
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physics responsible for the extra particles. 

In order to optimize the trileptons events, one should reduce the W ± Z background 
and ti background with suitable cuts. Those have been discussed in works [13J HU H5] , 
which focus on trileptons events from chargino-neutralino production, and after the cuts 
they find a total SM background of 19.6 fb. Here cuts designed in Ref [15] are as follows: 

1. 3 leptons with p T > 10 GeV. 

2. At least one Opposite Sign Same Flavor (OSSF) dilepton with Mqssf > 20 GeV to 
suppress low- mass 7*, J/^f, T, and conversion backgrounds. 

3. Lepton track isolation: p^trk < 1 GeV for muon and < 2 GeV for electron, where 
PtIa i s the maximum of any additional track within a cone R = 0.2 around the 
lepton. 

4. Ef ss > 30 GeV. 

5. No OSSF dilepton has invariant mass in the range 81.2 GeV < Mossf < 102.2 GeV. 

6. No jet with p T > 20 GeV. 

In addition to the SM background, we also comment on the background coming from 
ordinary superpartners. The GMSB model with stau NLSP may produce the trilepton 
events, e.g., xfxl ~~ tv t + ff — > Gtv t + Grf — >■ E™ 1SS + trilepton, and so on. Those 
MSSM backgrounds include tiny branching ratio BR(3r — > I + OSSF dilepton + E™ 1SS ), 
and hence the MSSM backgrounds could be reduced enough. 

Red dash-dotted and purple dotted lines in Fig. [3] show the cross sections for the 
direct production pp — > H'^H' + and pp —> H'^H'~ , respectively. They are calculated with 
the CalcHEP [16] implementing the Lagrangian Eq. ([1]) and the CTEQ6L code [17] as 
a parton distribution function. In the calculation, we took m^/o = m^,±, and set the 
center of mass energy to be yfs = 14 TeV. From the numerical result, production cross 
section can be parametrized as a(pp — > H'^H ) = 2.47 x 10 _4 (TeV/m^) 4 ' 1 pb, where 
m^/o = mjj,±(= mjj,). 

§ince their decay modes into two bodies are kinematically inaccessible, they decay into 
dileptons and H'l via off-shell weak gauge boson. Their branching ratios into dileptons 
are, therefore, uniquely determined; BR(ff$ -»■ llH'l) ~ 6.73% and BR(#' ± HvR'l) ~ 
21.32%. Here results have summed over electron and muon. Thus the 5a discovery reach 
for the trilepton signals is estimated as follows, 

Signal = cr(pp ->■ H§H f± ) x £ x BR(H'£ -> UH'£) x BR(if /± ->■ I'vH'l) 

\mjjj VlOOfb- 1 ; V 

Here £ stands for an integrated luminosity, and N signal (iV BG ) stands for the number of 
multilepton events (SM background events). Thus, assuming an integrated luminosity 100 
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fb -1 and demanding the SM background 19.6 fb , the superpartners of extra doublet for 
mjj, = mfj,± < 205 GeV would be discovered at the 5cr level. Indeed, for the discussion of 
discovery reach, we should mention the acceptance of detector. However it needs precise 
simulation of signal events, and is beyond the scope of this work. 

Next, we discuss the case that the lightest extra particle is singlet-like boson, S'. They 
are mainly produced by the decay of heavier doublet-like bosons accompanying the SM- 
like Higgs boson, h°. Accordingly, the signal events contain bb and the missing energy, 
and hence provide a distinguishable signature from the ordinary GMSB ones. 

The dominant background to the two 6-jets plus large missing energy events presum- 
ably comes from ti production. It can be reduced by the following cuts: 

1. Ef ss > 100 GeV. 

2. 6-jets with p T > 50 GeV. 

3. Ef ss + J2 E Tj > 1500 GeV. 

Here ^Tj indicates the transverse energy sum over untagged jets [38] ■ The most promis- 
ing event for the discrimination, therefore, would be missing energy events accompanying 
bb and energetic jet, pp ->■ HgH* -> h°S' + W ± S' -> bbS'+]etS'. In both GMSB models 
with stau NLSP and Bino-like neutralino NLSP, there would exists no background events 
of ordinary superpartners under the cut conditions. 

Gray dashed and red solid lines in Fig. [8J show the cross sections for the direct pro- 
duction pp — > H'°H' + and pp — > H'°H'~, respectively. The production cross section is 
parametrized as a(pp — > if /0 if /=b ) = 1.027 x 10~ 4 (TeV/rnf/') 4 ' 2 pb, where m#/o = m#'±(= 
run')- Since the branching ratio for H'° — > bbS' has a complicated dependency on model 
parameters, we take it as a free parameter. When charged Higgs boson is much heavier 
than W ± , H f± dominantly decays into W ± and 5", and hence the branching ratio of H ;± 
into jets and S' is BR(if /=t — > jets + 5") ~ 67.60 %. Thus 5cr discovery reach for them is 
estimated as follows, 

^signal = *(PP ->• H^H^) X C 

x BR(#'° -± bbS') x BR(iJ ,± -»■ jet + S') x (6-tag efficiency) 

~ 174 ( TeV Y' 2 ( C \ ( BR ( H '° ^ b bS')\ A-tag efficiency \ (37) 

WW V100 fb -1 / V 10% / V 50% x 50% ) 
>5 v / iV B ^. 

The cross section of the ti background for these events is given by [31], a t i = 0.89 fb 
(a t t = 0.72 fb) for the 2% (1%) b mistagging probability. Thus assuming C = 100fb _1 , 
BR(H'° —7- bbS') = 10%, 6-tag efficiency = 50% x 50%, and demanding the ti background 
1 fb conservatively, the extra doublets for m H m = m H i± < 260 GeV would be discovered 
at the 5cr level. 

Finally, we discuss the collider signature of the color-triplet states, introduced in order 
to maintain the successful gauge coupling unification, which is another key ingredient for 
the selection at the LHC. 
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Those colored extra particles are produced, and then they hadronize in the detector 
materials. The hadronized particles would be electrically either neutral or charged, and 
they can reverse the sign of its charge through the scattering in the detector materials 
[4"9| |50| I5T] . If a hadronized particle is electrically neutral and would not undergo the 
charge reversal, it leaves detectors. The resultant large missing energy without energetic 
tau lepton or photon would be the mark for the selection from ordinary GMSB model. On 
the other hand, some of the hadronized particles are charged, and they lose their kinetic 
energy through ionization with the detector materials. The ionization energy loss is a 
function of /?7 and the electric charge of penetrating particle [52]. When the energy loss 
and momentum of penetrating particle is measured, can be obtained, and hence its 
mass is determined. In addition, since massive long-lived charged particles produce a track 
in detectors, we can speculate its production rate. The estimated discovery potential with 
this method, however, is necessarily dependent on the scattering model, which of charge 
reversal predictions [49]. Thus the discussion of feasibility for the discovery requires model 
dependent full analysis, and we leave it for future work. 

5 Conclusions 

In this paper we proposed an extension of the GMSB model by adding extra doublets and 
singlet with parity odd under an additional Z 2 parity, while ordinary MSSM fields are 
parity-even. In this class of model, a natural candidate for DM appears as the lightest 
linear combination of additional Z 2 -odd fields. It has sizable interactions with nucleons 
and will be detected future /on-going direct detection experiments. We have also discussed 
typical collider signatures of this model and found that they may be discovered at 5 sigma 
level with an integrated luminosity of 100 fb -1 , depending on model parameters. 

The calculation of the relic abundance of bosonic-singlet DM indicates DM mass is 
around 250 GeV. If this is realized, this model could be confirmed by both the LHC 
experiment and future DM direct detection experiments. Furthermore, these DM can 
reproduce the DM-like events at CDMS-II. On the other hand, fermionic-doublet DM 
mass is around 1 TeV as shown in Fig. 1. In this case, it is very hard to discover these 
extra particles at the LHC, and DM direct detection experiments would not observe its 
signal. 
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